Intermittency is a consequence of turbulent transport in nonlinear systems 
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Turbulent flows transferring energy from a stirring 
range at large scales to the dissipation range at small 
scales consist of dissimilar components, a broad spectrum 
of eddies and randomly occurring, intermittent coherent 
structures. The first cascade leads to Kolmogorov-like 
finite-flux power spectra. The second component is par- 
ticularly visible in the anomalous short-scale behavior of 
the higher order moments of velocity differences. Com- 
pared to fully developed three dimensional turbulence, it 
is easier to quantify these two components in turbulent 
systems of weakly coupled dispersive waves. In this case, 
the Kolmogorov-Zakharov spectrum is a stationary finite 
flux solution of kinetic equations that follow from three 
and four wave resonances of the weak interactions. The 
second component emerges since the wave turbulence ap- 
proximation 1] is almost never valid at very low and very 
high wavenumbers where the 'weak' coupling approxima- 
tion breaks down, leading to the emergence of fully non- 
linear structures Q. In short, despite the fact that the 
amplitudes are, on the average, small, the weakly nonlin- 
ear dynamics can lead to intermittent and localized high 
amplitude events and anomalies in high order moments. 
It may also lead to a contamination of low order moments 
and to power spectra which, at least in some wavenum- 
ber ranges, are dominated by strongly nonlinear events. 
Nowhere is this more evident than in the illuminating 
studies of Cai, Majda, McLaughlin, and Tabak (CMMT) 
(later confirmed by Zakharov, et al. Q) on damped, 
driven and on freely decaying weakly nonlinear dispersive 
one dimensional wave systems. Indeed, CMMT found 
that, in damped and driven systems, there were some 
situations in which the nonlinear solutions dominated at 
almost all scales. In the freely decaying state, they found 
the Kolmogorov-Zakharov spectra were much more likely 
to appear. This of course was disappointing because the 
strength of the Kolmogorov-Zakharov solution of the un- 
damped, undrivcn kinetic equation is that it describes 
exactly what would be expected if an energy source at 
k = feeds at a constant flux rate through an inertial 
range to a viscous sink at k = oo. 

The purpose of this Rapid Communication is to demon- 



strate in a simple but representative model that in driven, 
damped systems in which there are fluxes of two con- 
served densities (energy and particle number), the real- 
ization of a statistically steady state demands the exis- 
tence of high amplitude coherent structures because they 
play a crucial role in balancing the energy flux budget. 
We show that a single cascade of weakly interacting waves 
transporting particles from a source to a sink leads to a 
steady net loss of energy of the system. Any steady so- 
lution requires contributions from the nonlinear terms in 
order to offset this energy loss. This is achieved by the in- 
termittent formation and destruction of high-amplitude 
structures. Fig.l shows the intermittent emergence of 
peaks for a discrete nonlinear Schrodinger equation 

iK + 0n+i + K-i - 20„ + 4>lK = H<t>, *) - V{<t>, t) (1) 

for a chain of complex oscillators. I? is a short- wavelength 
damping term while T drives the system on long space 
scales. Fig. 1(a) displays the typical long-time behavior 
of \4> n \ 2 in a sector of 30 lattice sites from a chain of 
N = 512 oscillators with periodic boundary conditions 
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FIG. 1: Numerical integration of the DNLS: (a) \(/> n \ 2 for a 
sector of 30 oscillators over 30 time units with strong damping, 
(b) Density of sites with the amplitude \(f>\ as a function of 
|</>| for weak and for strong damping and a Gaussian fit of the 
density for weak damping. 
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over 30 time units. \<t> n \ 2 is small on average, but high- 
amplitude structures emerge locally. The peak shows a 
'breathing' behavior, decreasing and increasing irregu- 
larly. Fig. 1(b) shows the average density p of oscillators 
with an amplitude \cf>\ for weak and for strong damp- 
ing forces. The density for small \<j>\ is Gaussian, but for 
larger \<f>\ the intermittent high-amplitude structures lead 
to a slower non-Gaussian decay. In the weakly damped 
system, the density increases slightly for high amplitudes 
so that there is a small hump near \(f>\ = 2. The number 
of these intermittent events increases with the driving 
force, but the hump remains in the vicinity of \<j>\ = 2. 
In the simulations we use time-periodic (S-kicks as damp- 
ing and driving forces T> and !F, which allows a simple 
control of the energy input and output: The driving in- 
creases the homogeneous mode cq — > (1 + Ajf)co, and 
the damping decreases modes on a short space scale as 
Ck — ► (1 - Xv)ck where c k = (l/N)J2 n 4'nexp(ikn) and 
the wavenumber k is in the Brillouin-zone [0, 2tt]. In 
numerical studies we apply the damping to the short- 
wavelength modes with 7ir/8 < k < 9ir/8, so that the 
impact of the damping does not decrease with the sys- 
tem's size. Analytical studies are simplified by restricting 
the damping to the mode k — tt. T> and T are zero in 
the intervals of one time unit between the time-periodic 
kicks. No important changes are found for shorter in- 
tervals between the kicks, or for continuous driving and 
damping. 

The dynamics between the kicks is governed by the 
nonlinear Schrodinger equation that derives as i(j> n — 
dH/difin from the Hamiltonian 

H = H2 -t- H 4 

= En 2 ^«^ n - </>n</>„+l ~ <A„<^n+l (2) 
- IV d> 2 (h* 2 

The modulus square norm or particle number 

A = Y,4>nK = Nj2c k c* k (3) 

n k 

is a second conserved quantity. In recent studies Q 
we have shown that this isolated Hamiltonian system 
with no damping and driving force forms localized high- 
amplitude structures as a statistical consequence from 
the thermalization under the constraint of its two con- 
served quantities. Such peaks are generated in a self- 
focusing process of low-amplitude waves with long wave- 
lengths and a low energy. A typical initial condition from 
which high-amplitude peaks emerge is the Raylcigh-Jeans 
distribution of the power \ck\ 2 = ((3(ui(k) + 7)) _1 with 
a positive temperature 7 is the chemical poten- 

tial, Lu(k) = 2 — 2cosfc is the frequency. As the system 
approaches its state of maximum entropy, the spectrum 
of low-amplitude fluctuations becomes flat so that the 
power is equipartitioned on all modes (p — > 0,7 — > 00). 
During this transformation of the low-energy spectrum 
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FIG. 2: (a): Spectrum for the weakly damped nonequi- 
librium system compared to a similar low-energy thermody- 
namic spectrum (f3(u>(k) + 7)) -1 of the corresponding un- 
damped and undriven Hamiltonian system, (b): Logarithmic 
plot of the spectrum for strong damping. The steps of the 
spectra at k = 7ir/8 result from the damping. 

to an equipartitioned spectrum, stable high-amplitude 
peaks with a negative energy emerge as a by-product of 
the production of entropy in the low-ampitude waves. No 
such peaks emerge from low-amplitude short waves with 
a high energy. 

In contrast, the long time behavior of solutions of the 
damped and driven system is governed, not by the values 
of H and A, but by the fluxes of both quantities. Each 
driving step T changes the amplitudes 4> n by A<j) n = 
AjfCo, so that the total number of particles increases 
by AjrA = £„<^A< + A<j> n <t>* n = 2A^iV|c | 2 while 
the damping changes the particle number by A-pA = 
— 2ApA r |c 7r | 2 for small values of \fj>- Consequently, 
there is a flow of particles from the source at k = to 
the sink at k = 7T, and the power spectrum (Fig. 2(a)) 
decays with the wavenumber. With most power gath- 
ered at low wavenumbers, the spectrum is similar to a 
thermodynamic spectrum (Fig. 2(a)) of the correspond- 
ing undamped and undriven Hamiltonian system with a 
low energy. The high-amplitude structures are generated 
just by the same mechanism in such an isolated Hamilto- 
nian system and in the damped and driven system where 
the permanent particle flow maintains the bias of the 
spectrum. The spectrum decays even exponentially for 
strong damping forces Fig. 2(b). Particle loss and gain 
are balanced when the flow in a sufficiently large system 
is constant so that 

AjrA + A V A = Q (4) 

or Ajf|co| 2 = Aulcjrl 2 . The driving parameter A^- and the 
damping parameter Ap regulate the particle number and 
the particle flux of the system. By choosing the driving 
parameter as Ajf = 2 n\c \' :! we orj tain a constant flow of 
particles A^A = AA into the system in the numerical 
simulations. The damping parameter Ax> is fixed. For 
a strong damping to forcing ratio, A-p = e~ 2 Ajr with 
e 2 <C 1, it follows from @ that \c^ \ — e|co| and 

A jc- 1 c 1 = eX-vlc-xl (5) 
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FIG. 3: (a): Change of H2 due to the damping and H4, due 
to the damping and driving force as a function of time, (b): 
Time average of the same energy flows as a function of the 
particle flow A A and the damping constant A -p. 

The input and output of particles also changes the 
quadratic coupling energy and the quartic energy. Gains 
and losses again have to match in a stationary nonequi- 
librium state so that 

Aj=H 2 + AjrHi + A V H 2 + A v Hi = (6) 

In order to understand the role of the strongly nonlinear 
terms in the energy flow, we analyze how it divides up 
among these four terms. The change of the quadratic 
coupling energy H 2 = N^Lo(k)ckC* k is given by the flux 
of particles times the frequency u>(k). The influx of par- 
ticles through the driving force leads to a zero energy 
influx AjrH 2 = w(0)A4 = since w(0) = 0. The damp- 
ing leads to a loss of coupling energy ApH 2 = — uo{-k)AA 
with Ct>(7r) = 4. Consequently, the particle flux leads to 
a net loss of coupling energy AH 2 = AjrH 2 + Ax>H 2 = 
— 4AA. This energy loss must be balanced by an en- 
ergy gain in the nonlinear term H4 for maintaining a 
stationary nonequilibrium state. The viscosity changes 
the nonlinear energy H A = — (1/2) ^n^ 2 as 

A V H 4 = 2X v NRe(d„c* 7r ) 

where d k — (l/iV)lZn \(j) n \ 2 4> n e lkn . The driving changes 
this energy contribution by 

AfHi = -2\rNRe{d c* ) 

Fig. 3(a) shows A?>H 2 , A^H^ and A^H^ as functions 
of time for a chain of 512 oscillators for a total particle 
flux per time unit A A = 0.512 and a strong damping 
to forcing ratio where A-p = 0.2. This corresponds to 
the simulation of Fig. 1(a), the strong damping case of 
Fig.l(b), Fig.2(b) and to Figs.4(c),(d). The fluctuations 
indicate single collapse events which have a smaller rel- 
ative strength for larger system sizes. Fig. 3(b) shows 
the time average of At>H 2 , A^H^, ApH^ as functions 
of the damping coefficient Xv and the particle flux AA, 
where the weakly driven case (AA = 0.512) corresponds 
to all other simulations. The impact of the driving force 
on the energy AfHn ss is negligible compared to the 
two effects of the damping AdH 2 and ApH^. This pre- 
dominance of Ax>H4 over A^H^ (i.e. XfReid^c*^) <C 



Ax>i?e(c? 7r c* ) ) follows from Ajr|co| <C A-d|cv| by equation 
(0 and from \d^\ ~ O(|do|)- The second relation is due to 
the spiky shape of \(j) n \ 2 (f) n that leads to a flat power spec- 
trum. In addition, the phases of d^ and c„ are correlated 
by the phase velocity of the peaks while doc^ oscillates 
randomly as it is governed by low-amplitude fluctuations. 
AjrHi yields a significant negative contribution only for 
high particle fluxes AA and weak damping forces Ap (at 
Xv = 0.02, AA = 2.5 in Fig.3b). 

For all other cases, both AjrH 2 and A^H 4 are zero or 
close to zero, so that @) reduces to Az>H 2 « —ApH^. 
This yields i?e(c*d 7r ) ~ A^c^c^ or the inequality 

Kl > 4|c T | (7) 

The importance of the nonlincarity A-pH^ in the en- 
ergy balance is reflected in the system's first three mo- 
ments: Defining = </>«(— 1)", the inequality J7J) reads 
lEK'nlVnl > 4|£Vn|- The first moment |£>n| ~ 
V AA increases with the particle flux. The second mo- 
ment IV'nl 2 = A is the particle number, which is small 
if the system is strongly damped and weakly driven. 
Equation (0) shows that the third moment exceeds the 
first moment, so that the distribution is strongly non- 
Gaussian. High-amplitude excitations (Fig. 1(a)) lead- 
ing to the non-Gaussian tails of Fig. 1(b) are necessary 
to maintain a stationary nonequilibrium state under the 
flow of energy and particles. 

To understand how the high-amplitude structures bal- 
ance the energy flow, we compute the input and out- 
put of energy and particles as functions of the amplitude 
1 0| of the oscillators. The driving kick feeds a num- 

ber AM\4>\)d\4>\ = En#J<|W<M+<iW 2ife(0„A^*) 
of particles to those lattice sites where where the 
amplitude is between \tp\ and \<j>\ + d\(j>\. The 
total particle gain as used in 0} is AjrA — 
J °° Ajra(|</)|)c?|(/)|. Similarly, the forcing changes the non- 
linear energy at these lattice sites by AjF/i4(|0|)<i|</!>| = 
En^i^i^i^i^i+dl^l 2 !^! 2 ^^*^^)- The increments 
of the damping Ax,a(\4>\) and Ap/i4(|</)|) are defined anal- 
ogously. Ajrh 2 (\(j)\) is again zero, and A-vh 2 (\<j)\) — 
Ax>a(\<t>\). Fig. 4 shows the time average of these 
influxes and outfluxes of particles and energy as func- 
tions of \(f)\. In the weakly damped case (Fig. 4(a), (b) 
with X-p = 0.02, AA = 0.512), the particles are fed into 
the system mainly at lattice sites with moderate ampli- 
tudes |0| « 0.5 (Fig. 4(a)), while particles are removed 
mainly at sites with high amplitudes \<fi\ « 2. While the 
input of particles has a negligible effect on the balance of 
nonlinear energy, the output of particles leads to a gain 
of energy peaked at \4>\ « 2. Remarkably, a very small 
number of sites with high amplitudes is responsible for 
practically the total outflux of particles and the corre- 
sponding feedback of nonlinear energy. 
This energy gain can be seen for the most simple case of 
an isolated real peak <j) n — x$ni whose height decreases by 
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FIG. 4: Average loss and gain of particles (a),(c) and of non- 
linear energy (b),(d) as functions of the amplitude \cj>\. (a),(b) 
is the weak damping case, (c),(d) is the strong damping case. 

Ax>(f> = Xt>Ctt <C x during the damping kick, so that the 
coupling energy changes by Api^ = — 8xAp0. This in- 
creases the nonlinear energy by Ax>H4 = 2x 3 Av4>- These 
fluxes are balanced if the peak amplitude is \ — 2 which 
agrees with the results of Figs. 4(a), (b). Particles removed 
from the tip of the peak have to be replaced constantly 
through the driving mode at k = 0. During the Hamil- 
tonian evolution, the focusing process increases the peak 
again while short-wave fluctuations are radiated. 
The input of particles in the strongly damped equation 
(Fig.4(a),(b) with \ v = 0.2, AA = 0.512) is maximal 
for low amplitudes \<f>\ ~ 0.4 similarly to the weakly 
damped case of Fig. 4(a). The damping leads to a sig- 
nificant loss of particles at \<j>\ « 1.5, and to a tiny loss 
of particles in the domain of low-amplitude fluctuations 
at \<p\ w 0.2, but, paradoxically, also to a particle gain 
at \cj)\ ps 0.7. While there is still a net loss of particles, 
this partial recycling of particles increases the ratio of 
energy gain per particle loss. This allows the system to 
satisfy the balance condition Q with a peak amplitude 
less than \<p\ =2. As a simple model for this mecha- 
nism we assume two peaks at the sites I and I + 3 where 
4>i = X > i — 4>i+3 are rea l an d positive. A damping kick 
decreases the higher peak to 4>i = x — At>4> and increases 
the lower peak to 0; +3 = £ + Ax></>. The decrease of the 
higher peak <\>i models the particle loss at \<$>\ w 1.5 in 
Fig. 4(c) and increase of the lower peak </>; +3 represents 
the gain of particles at |0| sa 0.7. The coupling energy 
loss is AdB.2 = —8(x ~ £) Ap(/> and the nonlinear energy 
gain is A23H4 = 2(x 3 — £ 3 )Ax></>. The gain of nonlin- 
ear energy can outweigh the loss of coupling energy for 
X>2/V3. 

In summary, coherent structures are an essential compo- 
nent of the transport of particles from the large scales 
to small dissipation scales. The particle flow leads to a 
steady loss of coupling energy, so that the fluctuations 



have a low ratio of energy per particle. Such long-wave 
fluctuations form high-amplitude structures and it is 
the role of the coherent structures to balance the energy 
input and output by exploiting the nonlinear component 
of the energy. To offset this energy loss, high amplitude 
structures containing a significant amount of negative 
nonlinear energy must be formed and destroyed. Energy 
is transfered from the nonlinearity Hi to the coupling H2 
during the peak growth, and pruning or destroying the 
peaks increases the systems total energy. The peaks are 
essentially the pipes through which energy and particles 
flow under the constraints of the balance conditions. 
Obviously, this effect follows purely from the dispersion 
and the nonlinearity, and it is not restricted to the one 
dimensional discrete system of our simulations. A similar 
behavior was found [2j for the two-dimensional continu- 
ous focusing nonlinear Schrodinger equation when the 
forcing is applied at a midscale wavenumber and the 
damping at large wavenumbers. This leads to an en- 
ergy flux to high wavenumbers, and to an inverse par- 
ticle flux to low wavenumbers. The particle flux builds 
condensates from which collapsing solitons emerge that 
carry particles to the dissipation scale and feed both par- 
ticles and energy back into the wave field at these high 
wavenumbers. With an almost Gaussian statistics of the 
waves and a Poissonian distribution of the intermittent 
coherent events, the system resembles a two species gas. 
The presence of collapses does not appear to affect low 
moment statistics of two- or three-dimensional systems 
if the damping is sufficiently strong, while the collapses 
contaminate the power spectrum generated by wave- wave 
interactions in nonintegrable one dimensional systems. 
The important point to stress is that no particles would 
be dissipated and no steady state would be achieved with- 
out strongly nonlinear collapses. 
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